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The quantum states filtering, whose general theorem was given by Bergou et al. (Phys.Rev.A 71, 
042314(2005)), should find it's important applications in present scheme, where we are trying to 
show that the problem of quantum states unambiguous discrimination may be solved by applying 
the argument of filtering. Let's use the quantum filtering, as an example, to show the basic idea of 
present scheme. Suppose there are N linearly independent states, if we are able to find a (N+l)- 
dimensional unitary transformation, R{u>) ( with uj is an adjustable variable(s)), which will be 
performed on each * > in the way like: R(w)\^!i >= dn{uj)\Dx > +f 1 (cu)\F >, R(u)\^ k >= 
Yll=2 dkl(w)\Di > +fk(uj)\F >, then, according to the definition of the operators for filtering, 
there should be: E^lo) = R' 1 ^)^ >< Dx\R{ui), E 2 (ui) = £)f =a fl^MI-D* >< D k\R(u) and 
E (lo) = R~ 1 (uj)\F >< F\R(u). With this this {Em{uj)} in hands, we could find the optimal 
operators which lets the function F(u)) = Y2i=i Vj < ^j|J^o(w)| 1 i r j >, with r\j to be the a priori 
probability of \<S>j >, have it's minimum value. For the system with N=3, there are three types of 
operations: (a) £?|* 2 >= £?|*3 >= 0, El |*i >= and E$ corresponds to fail; (b) >= 

if i ^ j for 1=1,2,3, and E b for failure; and (c) Ef|* a >= £i[#3 >= 0, £ 2 c |*i >= El\^ 3 >= 0, 
while E§ = I — 53t =1 Ek- We shall show that all these three types of operators, which may be 
performed on a N=3 systems, can be get by applying argument of filtering: the case a is in fact the 
filtering with N=3, case b can be viewed as successive filtering and the case c can also be solved by 
an argument of filtering in subspace. It can be shown that each case, which belongs to the above 
three, can be solved by reducing it to the problem of filtering. An important case of N=4 system, 
has also been discussed. 

PACS numbers: 03.67.Lx 



I. INTRODUCTION 



As a very recent development, the possibility of un- 
ambiguous discrimination between unknown quantum 
states can be potentially useful for many applications 
in quantum computing and quantum communications. 
The problem of unambiguous discriminating pure states, 
which are successfully identified with nonunit probabil- 
ity but witout error, was originally formulated and ana- 
lyzed by Ivanovic, Dicks and Peres [1-3] in 1987. Later, 
Jeager and Shimony solved the question of unambiguous 
discrimination of two known pure states with arbitrary 
a priori probability. Shortly after this result, Chcfics 
proved that only linearly independent pure states can be 
unambiguously discriminated [5]. The problem of dis- 
crimination among three nonorthogonal states was first 
considered by Peres and Terno [6], and the same ques- 
tion has also been discussed by Duan and Guo[7] and 
Sun et al. [8]. Chefles and Barnett also provided the op- 
timal failure probability and it's corresponding optimal 
measurement for a n symmetric states [9] , and an exper- 
imental set for discriminating four linearly independent 
nonorthogonal symmetric states was given by Jimenez 
et al. [10]. A new strategy for optimal unambiguous dis- 
crimination of quantum states was also offered by Ja- 
farizadeh et al. [11]. 

Unambiguous discrimination involving mixed state or 
a set of pure states, became an object of research re- 
cently. Several necessary and sufficient conditions for the 
optimum measurement have been given by Zhang et al. 
[12] and Eldar et al. [13]. Reduction theorems, which can 



simplify the discrimination theorem, have been developed 
by Raynal et al. [14-15]. Low bounds for the failure and 
the conditions for saturating the boumds, have also been 
studied [16-20]. There are only a few special cases have 
analytical solution for the quantum measurement, for ex- 
amples, the quantum state filtering [21-23], two mixtures 
with orthogonal or one-dimensional kernels [14-15], two 
mixtures in The Jordan basis [24] and other cases [26-30). 

In present work, we shall present a new scheme to 
solve the problem of quantum state unambiguous dis- 
crimination. Let's use the quantum state filtering orig- 
inated from [21-23], as an example, to show the ba- 
sic idea of present scheme. Suppose there are N lin- 
early independent states, the task of the quantum state 
filtering can be viewed as to find a set of operators 
{E m }, whose elements are defined by : E 2 \^i >= 0, 
Ei\^k >= with 2 < k < N, and E corresponds 
to fail. If we are able to find a (N+l)-dimcnsional 
unitary transformation, R(uj) with u> is an adjustable 
variable(s), which will be performed on each |^ > in 
the way like: R(u>)\^i >= dn{uj)\Di > +f 1 (u)\F >, 
R(iv)\* k >= £jI 2 dwMlA > +/ fc (w)|F >, then, 
according to the definition of the operators, there 
should be: Ei(u>) = ir^l-Di >< Di\R(u), 
= J2k=9 R~ l {u)\D k >< D k \R(w) and E (lo) = 
R^iw^F >< F\R(oj), With this {E m (oj)} in hands, we 
could find the optimal operators which lets the function 
F(u) = J2j=iVj < ^j\ E o(^)\^j >, with r)j to be the 
a priori probability of >, have it's minimal value. 

For the system with N=3, there are three types of oper- 
ations: (a) ^1*2 >= E?\V 3 >= 0, |*i >= and Eg 
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distinguishable outcomes of an operation is associated 
with linear transformation operators A mi 
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(2.1) 
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corresponds to 
introducing the 



with An leads to 
the discrimination 
reciprocal states >, which is defined as that which 
lies in TC, the N-dimensional Hilbert space for the N lin- 
early independent states \%Fj >, and is orthogonal to all 
for j ^ /, Chefles found that [5] 
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FIG. 1: An optical setting to realize the POVMs in an en- 
larged Hilbert space 



where \ej > form an orthonormal basis for Tt while Pj 
is the conditional probability, given that the system was 
prepared in the state >, that this state will be iden- 
tified, 



corresponds to fail; (b) Ef\^j >= if i ^ j for i=l,2,3, 



and E$ for failure; and (c) Ef\^ 2 >= Ef\^ 3 >= 0, 



>= ^il*3 >= 0, while Eq = I — El=iE k . We 
shall show that all these three types of operators, which 
may be performed on a N=3 systems, can be get by ap- 
plying argument of filtering: the case a is in fact the 
filtering with N=3, case b can be viewed as successive fil- 
tering and the case c can also be solved by an argument 
of filtering in subspacc. It looks as if each case, which 
belongs to the above three, can be solved by reducing it 
to the problem of filtering. An important case of N=4 
system, has also been discussed. 

Our present paper is organized as follows. Section II is 
a preliminary section in which we introduce the so-called 
double-triangle representation. In section III, we shall 
give a different way of solving the question of quantum 
states filtering. A concept of filtering in subspace will be 
introduced in Sec. IV. Two examples, discriminating three 
pure states and discriminating two mixtures for N=4, 
will be discussed in Section V and VI, respectively. In 
Sec. VII, we conclude the paper with a short summary. 



II. DOUBLE TRIANGLE REPRESENTATION 



A. preliminary 



Considering a quantum system prepared in one of N 
pure quantum states \&j >, where j=l, 2,..., N, if the 
states are non-orthogonal, no quantum operations can 
deterministically discriminate them. It is, however, pos- 
sible to device a strategy reveal the state with zero er- 
ror probability under the condition that these states are 
linearly independent [5] . Employing the Kraus represen- 
tation of quantum operations [31], each of the possible 



<TpM],A j ,\* j >=P j 6 jr . (2.3) 

In the terms of positive operator valued measures 
(POVMs) [31], the measurement can be expressed by 
defining the positive Hermitian operators 



Eh 



2i 1^4- >< tf-M 

|<VM->I 2 ' 3 jl ' 



e = A+A , 



(2.4) 



ith E () 



I, and it has been shown that 



the optimum measurement corresponds to the maximum 
eigenvalue of value of Ed = SjLi Ej being equal to 1 
[5]- 

Let G = >} for j=l, 2,..., N, the POVMs 

given above can viewed as one type of operations on G. 
There may be other types of operations on the same G, 
for example, if there are two known groups of states, 
d = >} for k=l, 2, K, and G 2 = >} 

for 1= L, L+l, N, G\ and G2 may have common ele- 
ments if L < K, we could also define a new set of POVMs 
{Ei, Eg} for i=l, 2, Ei can unambiguously tell whether a 
state |^? >e G belongs to d or not. Now, one may ask 
the question: could this {E m } also be expressed in terms 



of I*/- >< 



We shall give an answer to this ques- 



tion. According to the Neumark's theorem [32]: if each 
E m is an one dimensional positive operator, {E m } can 
always be realized by extending the Hilbert space to a 
larger space and performing orthogonal measurement in 
the larger space, while, as we shall shown, we are able to 
realize the discrimination of the quantum states just ac- 
cording to the definition of the operators, this fact makes 
it possible to read {E m } from their corresponding pro- 
jective operators in the enlarged space. We shall show 
how this basic idea works via the aid of Fig.l: the total 
space is defined to be H = H e + H A with { je^ >,\vj >}, 
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i=l, 2, N, and j=l, 2, M, for it's "in-space" while 
{\D k >,\Fi >} for it's "out-space" with 1 < k < K, 
1 < I < L, and N+M=K+L. H e is the Hilbert space 
where the states are defined: 



l*i >= 



N 
3 = 1 



e,- >, 



(2.5) 



H A is the subspace for ancillas, U(N+M) will couple this 
two subspace together. Let's use lo to denote the ad- 
justable parameter(s) in the unitary transformation, we 
can define \ip 
" out-space" 



>0 ut_ > and express it in the 



K 



fe=i 



with the normalization constraint Y^k=i \djk\ 2 + 
Yld=i \fji\ 2 = 1- K we want to unambiguously discrimi- 
nate all in G, we should find the general U{lo) which 
gives 



L 

H >° ut = da^Di > + >> ( 2 - 7 ) 

i=i 



then, ID, >< £),•! and 



Eli 



Fr >< FA should be 



the projective operator corresponds to Ej and Eq, re- 
spectively With the U(lo), we could write , for exam- 



pie, \Dj > in the "in-space" as \Dj > m = U 
\Dj > e +\Dj > A with \Dj > e and \Dj > A are two non- 
normalized vectors which lies in H e and H A , respectively. 
There should be 

L 

Ej{uj) = \Dj > e < Dj\, E (u>) = \ Fl >S< Fl \> ( 2 - 8 ) 



i=i 



and E Q (u) + J2f=i E j{ w ) 



With these operators 



in hands, we could get both the optimal operators and 
the maximum values for discriminating \^>j >. When 
the projective operators are expressed in the "in-space", 
there are written formally in terms of |e» >< ej\. If we 
could define {\ej >} by >} at the beginning, then 

we shall be able to complete the task of expressing E m 
in terms of >< ^f\- The argument above can also 
be generalized to other cases with different operations on 
G. 

In present works, we always adopt the proposals orig- 
inated from the works by Sun et al. [8]: (a) any pure 
state can be realized by a single-photon state and (b), 
according to Reek's theorem [34], any unitary transfor- 
mation matrix can also be realized by an optical network 
consisting of beam-splitters, phase-shifters, etc., all these 
optical elements should construct an one-photon interfer- 
ometer (OPI). The device in Fig. 2 is a typical four-port 
beam splitter which is used to realize a two-dimensional 
unitary transformation [^(w): 



— sin lo cos lo 
cos uj sin lo 



\u),r > 
\lo, d > 



\lo, I > 
\lo, u > 



(2.9) 



V t l)- 



CO 



\m,d) 



FIG. 2: A four-port beam-splitter used to realize a two- 
dimensional unitary transformation 

A property of this beam-splitter, which is frequently ap- 
plied in present works, should be noted: suppose there is 
an input 



(2.10) 



>= U 2 {uj)\<j> >= c Wtr \u> 
with the coefficients satisfy: 



>= c w j\lo,1 > +c WiU \oj,u >, 

parameters fo 
le output shou 

r > +c Wid \oj,d> (2.11) 



where c Wi i and c WiU are real parameters for simplicity 
after performing the ^(w), the output should be 



Cui, r — — Cu,,i sinw + c Wi „ cosw, 
Cui.d = Cu,i cosw + c Wi „sinw, 



if we choose 
sinw = 



coso> 



r 2 + r 2 ' 



(2.12) 
(2.13) 



(2.14) 



then \<j>' > = (c 2 j + c 2 u )\to,d > while it's output along 
\lo, r > is zero. 

B. which parameters are known? 

In present work, we shall deal with the case that all the 
states in G are linearly independent and their overlaps 
are also known. 

Definition 2.1: a N-dimensional matrix O(N) is de- 
fined by it's matrix elements 



(2.15) 



with constraint that \0%j\ < 1 holds for i ^ j. 

Certainly, O(N) is Hcrmitian. Using 0*(N) and 6(N) 
for it's conjugate matrix and transposed matrix, respec- 
tively, there should be O* = Oij = Oji. 
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Definition 2.2: A(N) is used to denote the adjoint ma- 
trix of O(N), Aij = a-ij, the inverse of O should be 



dct(O)' 



where det(O) denotes the determinate of O(N). 
Definition 2.3: tj is used to denote 

tj =< il>f\*j > . 



(2.16) 



(2.17) 



From the definition of the reciprocal states, if |*j > 
is a the reciprocal state of >, then exp{i</>j}\$! j > 
is also a reciprocal state of |*j >. We can always let 
tj =< ipj-^j >—< ipj\&f > by choosing a suitable set 
of {<j)j}. Defining 



I*" >= 



1 



i*j > -*.i* ± >). 



(2.18) 



one may verified that |*^- > and |*'! > form an orthonor- 



mal basis for |\I>j > and 



(2.19) 



it should be emphasized here that, either 1*+ > or |^>j >, 
is defined from all the states in G: 

Theorem 2.1: < *+|* fe >= tjS jk . (2.20) 

With the tj and O(N) defined above, we may intro- 
duced another transformation matrix: 

Theorem 2.2: denoting (<%) = I, and defining 



R(N) = 6(N)(-5 ij ), 



(2.21) 



there should be 

RiN)^ 1 - >= |* > . (2.22) 

Proof: formally, we can write |*j > as a linear combina- 
tion of {|*^ >} in the way like \^j >= Ylk=i h jk\^h >> 
there should be Oij —< *i|*j >= SfeLi ^jk < 

>= hjiU, which gives |*,- >= ELi^l*fc >• 
Let N=3, as an example, we have 

= |*1 > 
= 1*2 > 
= 1*3 > 

(2.23) 





+ 


< * 2 |* 2 L > 


+ 


< * 3 |* 3 L > 


o 12 \^i > 


+ 


l*^> 




032 |*3 > 






< * 2 |* 2 L > 


+ 


< *s|* 3 L > 


oial*] 1 > 


+ 


23 |* 2 L > 


+ 





<*l|* 1 L > < * 2 |*2 > < *3|*3 > 



It is possible to express j* -1 > in terms of |* > through 
introducing the inverse of R(N) 



R-\N) = ( tj 6. 



det(0) : 



(2.24) 



from Eq.(2.16) and Eq.(2.21) while the relation, det(O) = 
det(O), has been used [33]. Naturally, 

iT^AOl* >= |* X > • (2-25) 

Both R(N) and which are known from O(N), can 

be used to derive the value of tj. Let's use N=3, as an 
example, to give the derivation. From Eq.(2.25), we have 

ti 



det(0(3)) 
det(0(3)) 
det(0(3)) 



(aii|*i > +a 2 i|* 2 > +a3i|*3 >) = \^T>, 
(ai 2 |*i > +a 22 |* 2 > +&32|*3 >) = 1*2 >> 

(ai 3 |*l > +«2 3 |*2 > +«3 3 |*3 >) = |*3 >, 

(2.26) 



times < Vfe I on both sides of each equation, there are 

dct(0(3)) 

cralizcd to 



-i—ll — — 1 for j=l, 2, 3. This calculation can be gen- 



Theorem 2.3: tj = ^2™. (2.27) 

Now,we have shown how to get > from G, and 
their overlaps can be expressed thorough 

Theorem 2.4: defining the matrix O- 1 (N) by 



Olj(N) = oij=< *i-\*f>, 
there should be 



I ^1] 

ij 



(2.28) 



(2.29) 



Proof: we could suppose l^f 1 - > is known at first while 
|* > can be viewed as it's "reciprocal" state, and there 
should be 



6- L (i«y«)l*>=l*- L > 

to 



(2.30) 



by following the argument for the case where |^ > is 
known at first. Comparing it with Eq.(2.25), we find 



= (^ — (-^ (2.31) 



it can be written in the form of Eq.(2.29) by using 
Eq.(2.27). Some O^- shall be given in the appendix. 



C. the double-triangle representation 

A complete set of reciprocal states exists if, and only if, 
the state \ipj > are linearly independent while the recip- 
rocal states are also linearly independent, this fact will 



■5 



be used in deriving a set of normalized basis set {\ej >}. 
Letting 



1*1" > = |ei >, 



(2.32) 



1*1 > 



ei > + 



i2l 2 |e 2 >, 



the coefficients of I* 4/ > are decided by the two require- 
ments (a) it's overlap with > keeps unchanged and 
(b) the state should be normalized. These requirements 
may also used in deriving the coefficients of 1*3 >: sup- 
pose 



1*3 >=J24k\ek >, 



(2.33) 



fe=i 



from the requirements (a) and (b), there are three equa- 
tions 

<^|vl^> = 4, >=1, (2.34) 



their solutions should be 



l-\oi 2 \*4 2 , (2.35) 



c 31 = o 



13' 



°23 — 0r >i 



21"13 



L 32 



1 <S3 



ft 



(2.36) 



J- 12 



31 



C 32 I 2 " 



In principle, this process can be continued until we get 



all the coefficients, 



used as the matrix elements for 



the matrix C ± (N). Introducing another N-dimensional 
Matrix E, which is defined by E = (C -1 ) -1 , with it's 
matrix elements denoted by , we can define the basis, 
{\ej >}, in the way like 



|ei > - |*r>, 



e 2 > = 



K 2 P 



(2.37) 



After introducing this basis, every input state can be 
expressed in it by defining 



en =< e,-|*i > . 



(2.38) 



One may verify that there should be = if i > j 
according to theorem 2.1, this makes 



\e N >= \V N > 



(2.39) 



according to a simple reasoning, and the expression, 
\ e N >= Eili e ji\^i >> i s an equivalent form of it. The 



ft 
ft 



ft 
ft 



phase-shifter 



detector 



FIG. 3: the OPI for filtering 



matrix C with is a upper-triangle matrix while C 1 - is 
a lower-triangle matrix, for examples, 



cfi 

c 21 c 22 
c 31 c 32 c 33 



Cll C12 C13 

c 22 c 23 
c 33 




l*r> 

l*2 L > 





(2.40) 



(2.41) 



this is the reason why we call {\ej >} the double-triangle 
representation (DTR) . In the argument below, we always 
suppose that the states, either {\^>i >} or >}, have 

been expressed in the DTR. At the end of this section, 
we would like to emphasis again: if O(N) is known, then 
tj, ± (N), {|*+ >}, {|e 4 >}, C and C 1 - are also given 
at the same time. 



III. THE QUANTUM STATE FILTERING 

A. the POVMs for the filtering 

The quantum state filtering, which was termed in [21- 
23], is a special case of telling whether a state |*? > 
belongs to , G\ = {*i} or G 2 = {\^k >} for k > 2, 
with a non-zero probability of failure. The derivation of 
the optimal measurement strategy, in terms of {Ei,Eo} 
for i=l, 2, to distinguish |*i > from G 2 has been given 
and it is shown that this problem is equivalent to the 
discrimination of a pure state and an arbitrary mixed 
states. The quantum state filtering, as we shall shown, 
plays important roles in present works: (1) it's an excel- 
lent example to show how our scheme works while (2) 
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the filtering in a successive way will be used to complete 
other operations on G. The POVMs {E m } for filtering 
are defined by Ei\^ k >= for k > 2, E 2 \^i >= and 
Eq for failure, our task is to find an general unitary trans- 
formation U(uj) which transfers each state \&j > in the 
"out-space" in way like: 



(3.1) 



= dn|Di >+/i|F>, 

N 

|*fc >° ut = J2 dk ^ D i > +-^ F > ■ ( 3 - 2 ) 

This U(uj) can be realized by the OPI in Fig. 3. Keeping 
in mind that \^k > has no input along the rail e\, the 
input state should be \*S>\ > when the detector D\ has 
been triggered. By applying Eqs. (2.9-14), we are always 
possible to prevent the signals of |Wi > from appearing 
in the detector D k . Usually, a complex parameter, say, 
Cij may be expressed as 



\cij\exp{i(j)ij}, 



(3.3) 



with exp{i(f>ij} 



\dj\. In Fig. 3, a phase-shifter, $y, 



is placed in front of a beam-splitter denoted by u^, we 
always choose the phase-shifter 



while each beam-splitter takes the value 



sinwfei 



COS U)kl 





\cij\ 2 - 


c\ x sin 2 uj 




\cij\ 2 - 


c\ Y sin 2 w 




Clfe| 






ki, I 2 - 


2 • 2 

cf j sin uj 



Cll 



o 21 ti 



with wn = w and =< ej|* >, for examples, 
*i, C12 = 
0, c 22 = 



C21 



Vi-I^l 2 ' 

*2 



(3.4) 

(3.5) 
(3.6) 

(3.7) 
(3.8) 



VI-K2I 2 
In Fig. 3, we could read 

\D! > m = expj-i^n}^ 1 ^!!)!^! >, (3-9) 
and get 

\D 1 > in = -sinw|ei > + sinw|i; >, (3.10) 

where $n = and Eq.(2.9) have been used. Through a 
similar argument, we could arrive at 



COSW21 smuj\v > 



\D 2 > m = \D 2 > 
\D 2 > e = - sinw 2 ie 

it can be proved that 

-cos we* 2 



(3.11) 

e 2 > + COSW21 cosw|ei >, 



\D 2 > e 



1*9 >• 



cos 2 uj + sin 2 uj I | 2 



(3.12) 



In fact, we may use the relation, E 2 (uj) = I e — E\(oj) — 
E (uj), instead of giving all \D k > e in detail. Using 
Eq.(3.1) and Eq.(3.11), we may get 

|Di > e - -sin^l^j 1 >, (3.13) 
|F>e = (I e -|Di> e <Di|)|*i> 

V< *i|(I e -|Di > e <£>i|)|*i >' 

and the POVMs of filtering should be 

E 1 (uj) = !£>!>«< Di|, £ M = |F> e <F|, (3.14) 

E2(W) - Y.\ Dk > e <^l =I e -£lM-W 



fc=2 



If the POVMs {i^m} were known, then the calculation of 
the optimal value of filtering should be easily completed. 
Suppose rjj is the a priori probability of >, we denote 
Pif(N-i) ( w ) an d ^/(jVi) ( w ) the average value of success 
and failure of filtering, respectively, 



Pi/(jv_i)(w) - ^Tr{E lPl 

i=l 

^i/fjv-i)^) = Tr(£ p), 



(3.15) 



AT 



fi/(JV-i)M + F 1/(JV _i)M = ^Jfo, (3.16) 



with pi = 771 l^i >< Vi|, P2 = Y,k=2Vk\^k >< *fe| and 
p = pi + p 2 - A simple calculation shows that 

N 

F 1/{N _ 1) (u;)= m \f 1 \ 2 + J2m\fk\ 2 , (3.17) 

k=2 

with /j =< F\ipj > and 

(3.18) 



fi = V 1 ~ ^sin 2 / fc 



Olfc 



2 • 2 

z sin w 



The optimal value of i*i/(jv_i)(w), F°J^ n _ 1 ^(uj) = 

Fi/(N-i)( UJ ° pt ), is defined to be minimum value of 
Fi// n _i\(oj) in the domain of < uj < 2ir. From 
Eqs.(3.17-18), there is 



dF 1/{N _ 1) {uj) 2 . 

— t — = K sin 2cj(?7i 

aw 



Efc=2 %Kfcl 
1 - i 2 sin 2 w 



(3.19) 



and F°J* n _ 1 ^(uj) happens at dF 1 /( N _i)(uj) / duj\ 

0. Now, we are able to give the optimal values of filtering: 

( a ) if fMk\ 2 > 1, by letting sin 2 ^' = 0, we 

have 



JV 



F 1 7 ( V-i ) M = m + El°^l 



(3.20) 



fc=i 



7 



(b) if i - 1\ < ^Ef^ tl° lfe l 2 ^ x > throu s h lcttin s 

(3.21) 



sin a; 



X ^'o lfe | 2 



=2 m 



t 2 
'-I 



we arrive at 



^/(V-i)^) ~ 2 



\ 



N 



(3.22) 



fc=2 



and (c) if \jYl!k=2 ^fl !*! 2 < 1 — if, the optimal value 
should be 



EN I |2 

_ fc=2 Vk\Olk\ 

i/(jv-i)1> w J - '/H 1 " L i> "I" f—^2 > 



(3.23) 



while sin 2 w P* = 1. Substituting w°p* for w in Eqs.(3.13- 
14), we could also get the optimal POVMs 




FIG. 4: the OPI for filtering with N=2 



B. an example: filtering for N=3 

The filtering of N=3 is a case with a fully analytical 
solution and an optical implementation of the optimal 
strategy [23], we shall show, via a simple optical setting, 
how to recover all the optimal values in [23]. For N=3, 
there is 



1-t 



Defining 



2 _ |oi2 | 2 + |oi3 1 2 - 012023031 - O21O13O32 



1 - I O23 1 



A = ry 2 |oi 2 | 2 + %|oi 3 | 2 , 



. (3.24) 



(3.25) 



we can write the general results of filtering to the N=3: 
(a) for A > rfi, there is 



F$ = m + A, 



(b) for ?7i(l - t\) 2 < A < r/i, there should be 

F<[% = 2^A, 
and (c) else, A < 771(1 - i 2 ) 2 , 



(3.26) 



(3.27) 



(3.28) 



Besides all this optimal results, we could also get the 
optimal POVMs for filtering with N=3. 



C. filtering with N=2: discriminating two pure 
states 

The problem, how to discriminate |*i > from |* 2 >i 
is the most interesting case in the field of quantum states 



unambiguously discrimination. Here, it can be solved as 
a special case of filtering with N=2. The present solu- 
tion is complete in the sense that: not only the optimal 
values but also the optimal POVMs should be given at 
the same time while the POVMs have the forms required 
by Eqs. (2.2-4). The OPI in Fig.4 is for the filtering with 
N=2 and it takes < *i|* 2 >=< *2|*i > for simplicity. 
From the DTR for N=2, the basis vectors should be 



|ei >= |*i >, |e 2 >= |*2 >, 
and the states have the forms 



(3.29) 



|*i > = V 1 - |oi2| 2 |ei > +o 2 i|e 2 >, 

1*2 > = |e 2 > ■ (3.30) 

With known parameters for N=2, which have been given 
in the Appendix, we have 



IVr > = 
Hi > = 



1*1 > -02l|*2 > 



V/1-|012| 2 
1*2 > -Ol2|*l > 

V 1 - I012I 2 



(3.31) 



by applying Eq.(2.25). The POVMs for discriminating 
two linearly independent states, |*i > and |*i >, should 
be 



E x = sin 2 oj\^>i >< ipt\ 



E-2 



Eq 



COS 2 U) 



(3.32) 



cos 2 to + sin 2 w|oi 2 | 2 
(|*i > -i 2 sin 2 w|* r >)(< *i| - t 2 sin 2 w < *^ 
1-t? sin 2 u> 



while tx = \Jl — |oi2 1 2 , and 



F V i(a;) = m (l- t{ sin 2 ui) + m ]° 12l 2 ■ 

1 — q sin ui 



(3.33) 
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Now, we could give the optimal values and the opti- 
mal POVMs at the same time: (1) for .f^ > -r^-,, let 

y ' \ VI |012 1 ' 

sin 2 uj opt = 0, there should be 



E opt = 0, E opt = >< 
E° pt = |*i><*i|, 

F l7l( W ) = Vl + V2\0l2\ 2 , 
P$(U,) = V2(l-\o 12 \ 2 ), 

(2) if |o 12 |<yf<^, by letting 

1 - 



(3.34) 



sin 2 u opt 



012 



l-\o 12 \ 2 
we shall get the optimal POVMs 

1 - 



(3.35) 



"72 



£ 2 opt 



E opt 



Mm 



1 



l°12| 

/ail 



\H ><H\ 



\H><H\, 



l-|o 12 | 2 



3 = 1 



which give the optimal results 

f$(u) = 2Vmm\oi2\, 

2 
3=1 

and (3) when . h 2 - < 1 012 1 , , through choosing sin 



(3.36) 



(3.37) 



1, we arrive at 

E opt = X tf r |, E opt = |* 2 X * 
while E?, pt = 0, all these operators leads to 



21, 



V 

jopt 



m\Ol2\ +V2, 



2 ,,,opt 



(3.38) 



(3.39) 



p° p l{u3) = m (i-\oi2n 



Compare with other methods of solving the same ques- 
tion, the present scheme states that the general POVMs 
may be given before the decision of the optimal values of 
success and failure. 



IV. FILTERING IN SUBSPACE 

A. filtering with the background 

Suppose Gi = >} and G 2 = {|*; >}, a new 

operation of G can be specified by the definition of 



CO 



FIG. 5: The OPI for discriminating two states with the back- 
ground 



the POVMs as #i|* 7 >= if |* 7 > belongs to G 2 , 
E 2 \^>/ >= if I*? > belongs to G\ and E$ corresponds 
to failure. If a state, say, \^f g >, is shared by both Gi and 
G 2 , then E\\^> g >= E 2 \^ g >= according to the defi- 
nition of the POVMs {Em}. We call this case the name 
of discriminating with the background. In this section, 
we shall consider a simple case of discriminating > 
from > with |* 3 > as the background. In the DTR, 
this operation on G can also be viewed as a filtering in a 
two-dimensional subspace. 

For N=3, the the basis vectors in the DTR should be: 



|ei > = > 
|e 2 > 
|e 3 > 



*2 > -°i 2 \^i > 



*3 >, 



^2? 



(4.1) 



with 



°12 = 



O13O32 — O12 



VC1 - |023| 2 )(1 - |013| 2 )' 

while the matrix C(3) takes the form 



(4.2) 



h °t: i3 °! 031 



V 1 - |o 32 | 2 



with t\, which holds for N=3, is 



032 
1 



(4.3) 



h = 



/det(0(3)) 



1 



1 023 1 



(4.4) 



The U(u>) in Fig. 5 is required to transform each \^S>j > 
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to |*j > out as 

|*i > out = di\Di > +f 11 \F 1 > +fi 2 \F 2 >, 

|* 2 > out = d 2 \D 2 > +f 2 i\F 1 > +.f 22 \F 2 >, (4.5) 

|*3 >° ut = \F 2 >, 

and this goal can be reached, as we shall shown later, 
by applying the argument of filtering. The J7~ 1 (w) shall 
give 

l-Di > e = -sin^l*] 1 >, 

\F 2 > e = |*3 >, (4.6) 
(7 e - \F 2 > e < F 2 \ - |d > e < DxDI*! > 



1*1 > e 



v/i-|/i 2 | 2 -Mil 5 



with |L> 2 > e has given by Eq.(3.12). The POVMs are 
defined by 

Ex{u) = |Di> e <Di|, E 2 (lo) = \D 2 > e < D 2 \, 

E {lo) = \F X > e < Fi| + |Fi > e < (4.7) 



Defining 



FrH = J>|/,i| 2 , 

3 

f 2 h = 5>|/ i2 | 2 , 



(4.8) 
(4.9) 



i=l 



with fij =< Fj\^i >, the average value of failure should 
be 



F(w)=F 1 (w)+F 2 (w). 
A simple calculation shows 

2 

F 2 {lo) = ^^|o 3i | 2 + r] 3 , 



(4.10) 



(4.11) 



it is still in a special form of filtering with N=3, see 
Eq.(3.26). The F^uj) is left to be decided by the gen- 
eral results of filtering with N=2. 

Formally, |*i > and 1*1 > can written by 



|*i > = V 1 - I c i3| 2 l*i > +ci 3 |e 3 >, 

1*2 > = y/l - |c 23 | 2 |* 2 > +c 23 |e 3 >, (4.12) 

with 1*^ > and |* 2 > are two normalized states defined 
in the subspace specified by {\ej >} for j=l,2: 



, , _ cnjei > +c 12 \e 2 > 

' X> " Vl-lcial 2 ' 
|*' 2 > = |e 2 >, 



and their overlap should be 
o'i2 =< *il*2 >= 



012 — Ol3°32 



v /(l-|03l| 2 )(l-|032| 2 ) 



(4.13) 



(4.14) 



Now, in the two-dimensional subspace with {|ei >, \ e 2 > 
}, our task is to discriminate |*! > from |* 2 > with 



v' 3 =V^-\o 3j \ 2 ), 



(4.15) 



for j= 1,2, to be their a priori probability, respectively. 
According to our discussion about filtering, we have 



= V [(l-(l-\o' 12 \ 2 )sm 2 uj) 

|2 



+ 



%l°12l 



l-(l-| ' 12 |2)sin 2 ^' 



(4.16) 



which is equivalent with the one given by Eq.(4.9). This 
is the reason why the present case is viewed as a process 
of filtering with N=2, certainly, in the subspace without 
the background. It's optimal results have nearly the same 

forms for filtering with N=2: (1) for 



sin 2 w opt = 0, there should be 



| > rco, let 



F opt {uj) 
P opt {uo) 



F 2 ( W ) + r ? i+r ?2 | ' 12 | 2 , (4.17) 



^ 2 (i-K 2 l 2 ), 

(2) if K 2 | < 7f < j^, by letting 

i _ . fu£ 



sin 2 co opt 



12 



*12 I 2 



1 - W 

we shall get the optimal results 

F opt (uj) = F 2 (w) + 2^f 2 \o' 12 \, 



(4.18) 



(4.19) 



3 = 1 



12 h 



and (3) when y < |°i 2 |> through choosing sin 2 uj opt 
1, we arrive at 



F°*{u) = F 2 (u;)+ V [\o[ 2 \ 2 + m , 

p° p \u) = »?i(i-K 2 | 2 ). 



(4.20) 



It should be noted that F 2 (uj), which has been given in 
Eq.(4.11), is in fact a constant. The present argument, 
which is suitable for discriminating G\ = {\if>i >, |* 3 >} 
and G 2 = {\ip 2 >, |* 3 >}, can be generalized to the 
discriminating two general mixtures sharing part of states 
in comm. 



B. discriminating two mixtures in Jordan basis 



Suppose there are two mixtures, 

K 



i=l 
2K 



P 2 = J2 Pjl*i><*il 
j=K+l 



(4.21) 
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H'l 



<5 



^5 



FIG. 6: The OPI for filtering the mixtures in Jordan basis 



with Yn=\Pi = Y^K+iPi = l - Lct to be the 
a priori probability for pk, k=l, 2, we may introduce 
i]i = \\Pi as the a priori probability for ["J^ > in 
G\ = >} while r^j = \2Pj as the a priori probability 
for > in G 2 = {\^j >}. If < * fe |*fe/ >= for k ^ k' 
except < ^il^i+fc >= cos^, pi and p 2 are called in Jor- 
dan basis. Defining the POVMs: E^Vj >= E 2 \^ l >= 
while Eo for failure, this {E m } can be get, as it has 
been show by the works in [21-23], through discriminat- 
ing pairs of pure states in each subspace. 

The OPI in Fig. 5 is used to discriminate these two 
mixtures, pi and p 2 . The total Hilbert space here is 
defined by H — J2iLi each H l is a two-dimensional 

subspace with it's basis as 



|*i > -COS^|*K+i > 



sin ( 



the two states in this H e are 



(4.22) 



|*i > = sin^|ei > + cos0j|ej, >, 
\«S> K+i > = \e 2 >, (4.23) 



while their reciprocal states 



1*+ > = 



\*K+i > = 



\Vi > -cos9i\y K+l > 
sin 0i 

\^K +l >-cos6 l \^ l > 
sin 9j 



(4.24) 



In this H l , our task is to filter \*S?i > from l^^+i >, the 



POVMs, {E l m }, to complete this task should be 



COS 2 U>i 



cos 2 oJi + sin un cos 2 6{ 

2 

fe=i 



(4.25) 



the average value of the failure in H l should be 
F^LOi) = 7^(1- sin 2 0; sin 2 w;) 

T]K+i cos 2 9i 



1 — sin 6j sin uj. 



2 ,.' 



(4.26) 



it's optical values are given by the theorem of filtering 
with N=2. Finally, we can define the POVMs {E m } by 



K 



(4.27) 



for m=0, 1,2. The average value of fail can be expressed 

by 



K 



i=i 

while it's optimal value 



(4.28) 



(4.29) 



i=l 



where oj° pt should depend on the actual value of the pa- 
rameters, rji, rjK+i and cos#j, here, this requirement has 
also been pointed by the recent work [20]. 



V. THE SUCCESSIVE FILTERING FOR 
DISCRIMINATION OF PURE STATES 

A. the optical realization of U(uj) 

In present section, we shall show the POVMs, which 
are defined for discriminating of pure states, can be real- 
ized in an enlarged Hilbert space by applying the succes- 
sive filtering. The OPI in Fig. 7 is designed to discrimi- 
nate three linearly independent states, \^>i > for i=l, 2, 
3, here. The U (oj) realized by this OPI can be written as 



U{w) = R{2)R{1), 



(5-1) 



with R(l) denotes the unitary transformation done by 
the beam-splitters and phase-shifters, Uji and $ji for 
j=l, 2, 3, on the left part of Fig. 7, while R(2) denotes 
the unitary transformation realizes by, uiki and <!>fci for 
k=2, 3, the beam-splitters and phase-shifters on the right 



11 



{H3 



and in this run of filtering, the average value of the failure 
should be 



F 2 (uj) = r,' 2 (l - (1 - |o' 23 | 2 ) sin 2 wa) 



+ - 



^|o' 23 | 2 



l-(l-| ' 2 3| 2 )sin 2 c 2 ' 

with U>22 = W 2 . If 

V2 = m = v, I/21HI/31 1, 

then we shall get 

F 2 pt = 2?7|o23 - /21/31I, 

by letting 

sin 2 w opt_Iz_K3i 



(5.8) 

(5.9) 
(5.10) 

(5.11) 



^23 I 



FIG. 7: the OPI for discriminating three pure stats 



B. the POVMs for discriminating three pure states 



part. At first, the R(l) is defined to filter |*i > from the 
states, |*2 > and |*3 >, here. 

> = dn|Di > +/n|Fi >, (5.2) 

> = Vl-l/fcil 2 !^ > +/ fc i|fi >, 

for k=2, 3, while 1*^ > lies in a two-dimensional sub- 
space H', which is specified by it's basis as {\u>ki,r >} 
for k=2, 3, and 



1*2 > = \A - l°2 3 l 2 l^2i, r > +</ 32 |*3 >. 
\% > = e^\u 3 i,r >, (5.3) 
where o' 23 =< * 2 |* 3 >> using Eq.(5.2), we get 

023 - /21/31 



Ono = 



(5.4) 



23 V(l-|/2il 2 )(l-|/3i| 2 )' 
In this run of filtering, the average value of fail should be 

3 

with fji =< Fi\^j >. After the first turn of filtering, 
we are left with two states, |* 2 > and > in H', with 
their a priori probabilities to be 



Vk=Vk(l-\fki\ 2 ), 



(5.6) 



for k=2, 3, respectively. R(2) is designed to filter \W 2 > 
from I \I>3 > in the way like: 



R(2)\% > 



?(d 22 \D 2 > +f 22 \F 2 >), 



VI-I/21I 2 

fl(2)|* 3 > = 1 =? (rf33|£>3 > +/32I-F2 >), 



V/I-I/Sll 2 



(5.7) 



By performing the R(2) after R(l), the state |*j > are 
transformed into: 



|*i > 
1*2 > 
1*3 > 



out 



out 



out 



du\Di > +fu\F 1 >, 

d 22 \D 2 > +hi\F! > +f 22 \F 2 >,(5.12) 

d 33 \D 3 > +/3l|^l > +f 3 2\F 2 > . 



With U-^w) = R- 1 (1)R- 1 (2), we can arrive at 
\Di > e = - sinwi|'0] L >, 



\D 2 > e 



1*1 > e = 



coswi sin lo 2 



,J cos w 2 + sin 2 u)i\o 
|*i > -ti sin 2 wl*] 1 > 



± 12 
12 1 



1*2^ >, 



(5.13) 



|F 2 > c 



'l — i 2 sin 2 wi 
(I e -|J 2 > e < £> 2 |-|*i > e < ^1)1*2 > 

Vl-|d 22 | 2 -|/2l| 2 ^ ^ 



with djj =< -Dj|*i >, /ij =< i^|*i > and U\ — Wn. 
For discriminating three pure states, the POVMs are de- 
fined by 



£7i(w) = |Di > e <£i|, 

£ 2 (w) = |L> 2 > e < D 2 |, 
2 

£oM = > e <F fe |, 



(5.14) 



fc=i 



£7 3 (w) = I e -^^-So 



1=1 



the average value of failure is defined by 
3 2 
F(u) = < >\ 2 + ^2m\< *i\F 2 > | 2 , 



(5.15) 
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it can be proved that 

F(u) = F 1 (u)+F 2 {u) 



(5.16) 



the F 2 (oj), which has been given in Eq.(5.8), is in the 
form of filtering with N=2. With calculations that 



hi = y 1 - qsinwi, 

, Olfc 

Jfei = 



\J\ — t\ sinwi 
for k=l,2, should be 



Fi(u) = 771(1 - tisinwi) + 



Efc=2 %|Olfe| 



(5.17) 



(5.18) 



1 — 1| sinwi 

certainly, it is also in a typical form of filtering with N=3. 



and (2), is si < 2s 2 , the optimal value 



I, si 



F opt = U— + 2s 2 ), 
A s 2 



with the choice of 



• 2 opt 

sin av — 



1 - ^ 

S2 



1 - 



2a| 

1 + S2 



(5.25) 



(5.26) 



substituting w° pt for Wi in Eq.(5.11), which is u>\- 
dependent, we can get the actual optimal setting for w 2 p *. 
One check that: the optimal values for F(lu), which have 
been given in Eqs.(5.24-25), are consistent with the op- 
timal results in [8]. 



VI. THE SUCCESSIVE FILTERING FOR 
DISCRIMINATION OF TWO MIXTURES 



C. the analytic optimal results for a special case 

Usually it is difficult for us to give an analytic solution 
for the optimal values of the F(oj) in Eq.(5.16), while 
the following case, which has been discussed in [8], is an 
exception. Considering the case, where < ^ 1 1 >=< 
^11^2 >= s i an d < *2|*3 >= s 2 under the conditions 
that s\ < s 2 , we find that 



/det(Q(3) 
I- si ' 



det(0(3)) = 1 - 2s\ - s\ + 2s\s 2 , (5.19) 

according to the results given in the appendix. Suppose 
i]j = 3 for j=l, 2, 3, we could find 



/21 — /31 — 



si 



y/l — t\ sinwi 
Using Eq.(5.10), there should be 

F 2 pt {ui) = ||s 2 - 



1 4-2 • 2 

1 — if sin uj\ 



(5.20) 



(5.21) 



the optimal value is defined as the minimum value of the 
function 



F(u) = F 1 (LJ 1 ) + F% pt (u; 1 ), 



(5.22) 



and it should depend on the actual situations about the 

2 

Oif (1) if 2si < 2i by letting 



• 2 opt 

sin uj{ 



1 - 2si 

o .0. 

1 - 



1 + S2 



we shall get the optimal result 



F o P t = _ (2si _ Sa)) 



(5.23) 



(5.24) 



A. the optical realization in the enlarged space. 



Suppose there are two mixtures 



Pi 



2 

£ 

i=l 



Pi\^i >< 



Pi 



3=3 



(6.1) 

with 5Z i=1 Pi = Ej=3 Pj = 1> each with it's a priori 

probability to be Ai, and Ei=i ^« = 1- Letting rji = \\-Pi 
for i=l,2, and i]j = X 2 ■ Pj for j=3, 4, the above question 
can also be viewed as an operation on G = {\^k >} with 
r/k as it probability for k=l, 2, 3, 4, we are required to tell 
if a state |* ? >G G belongs to G x = {|*i >, |* 2 >} or 
G 2 = {|^3 >, ^4 >} while there is a non- vanishing prob- 
ability for fail. In terms of POVMs, E^ipj >= E 2 \^i >= 
holds for \^fj ><E G 2 and |$, >€Gi, respectively. Cer- 
tainly, there is E = I e - ^Li E k- Thc O pl in Fi g-8 is 
designed to realize the unitary transformation in the way 
like 

U{u) = R(2)R(1), (6.2) 

where R(l) is the unitary transformation for filtering 
|5 , 2 > from the rest of the states in G, R(2) is used 
to filter |*i > from |* 3 > and |* 4 >, whose definitions 
shall be given later. It should be noted that, when the 
detector D\ fired, we can not tell whether this signal is 
from |*i > or \^f 2 > since the fact that these two states 
may have non-zero coefficients, C12 and c 22 , along the 
rail e 2 , respectively. It is certain that this signal can not 
come from the states, |*3 > and |\T/4 >, according to 
our discussion of DTR. One may compare the present 
R(l)\^j >, with the one in discriminating three pure 
states, 



i?(l)|* 2 > = d 2 l|£>l >+/ 2 l|Fl > 



l/i 



Mnl 2 l*i> 



+d 11 \D 1 > +/n|Fi >, 
v/l-|/ fe i| 2 |f' fe > +/ fc i|Fi >, 



(6.3) 
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for k=3,4. Defining a Hilbert space H' with it's basis as 
{1^21, t >, |o»3i,r >, \e\ >}, there are three states: 



l*i > = 



Cll 



v/l-|/n| 2 -|rfii| 



ei > 



+ 



y 3 1 "34 "41 
'1-K 3 | 2 



w 2 i,r > +04i|^4 > 



1*3 > = ^l-|<3l 2 l^i,r>+o 4 3l*4> ; (6-4) 
|<> = e^\ujsi,r>, 

with o-j =< >, these overlaps can be derived from 

Eq.(6.3) by the requirement that the keeps unchanged 
when the unitary transformation is performed on the in- 
put, 



°ik = 



o\k — /11/fei 



v/(i - I/11I 2 - Mnl 2 )(i - l/feil 2 )' 

for k=3,4, while 

034 — /31/41 

v /(l-|/3l| 2 )(l-|./4l| 2 r 



°34 = 



(6.5) 



(6.6) 



In principle, we could realize R(2) as the filtering for 

N=3, 



i?(2)|*i > = 



R(2)\*' k > = 



di 2 \Di > +.fi 2 \F 2 > 
^a-l/nP-ldiil 2 )' 
Et=3dki\Di > +fk2\F 2 > 



• (6-7) 



In the first run of filtering, the average value of failure is 
defined by 



^H = f>,|/ 31 | 2 , 

while the one for the second run of filtering is 
F 2 {lo) = 771(1 - c' n sin 2 w 2 ) 



(6.8) 



4 U' 1 2 



+ 



lfel 



1 — c' n sin 2 W2 ' 



in a standard form of filtering with N=3, and their 
a priori probability is 



r?i= ? ?i(l-|/ii| 2 -M 11 | 2 ), 
Vk = Vk(l - l/feil 2 ), 



(6.10) 



for k=3, 4, respectively. The coefficient, c' n , may be read 
from Eq.(6.4) as 



{H3 



OHZP-D 



FIG. 8: The OPI for discriminating two mixtures 



B. the POVMs realized by the OPI in Fig.8 

As we have shown, it's possible for us to get a U{lo) 
which transforms each > in the out-space: 

out 



*2 >° Ut = d2l|Dl > +/21IF1 >, 



|v&i> o Ut = Y^d lj \D j >+Y,hi\Fu>, (6.12) 

i=i i=l 

4 2 



Using these expressions and the inverse R 1 (l)R x (2), 
we can get 

| A > e = — sincji|e 2 >, 

IP (I e -|^l > C < A|)|*2> 

|Fl> " 7i = CT ' 

\D 2 > e = - sinw 2 |ei >, 



(6.9) \F 2 > e 



(6.13) 

(l e - Eli | A > e <A|-|Fi> e <Fi|)|^i> 

i-|/iiI 2 -ELiMhI 2 ) 



with dij =< Dj\$>i > and f tj =< F j \9 l >, the POVMs 
arc defined by 



E 1 = ]T|A > e < A|, 
2 

So = ]T|F, > e < Fi\, 



(6.14) 
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FIG. 9: the optimal results of F opt in Eq.(6.16), the present 
results are described by the solid line, while the analytical 
results in dash line. 



FIG. 10: the optimal choice of ui 



two mixtures: 



With these operators in hands, we could define the aver- 
age value of fail 

JV=4 

F(u)='£i lj <* j \E \* j >, (6.15) 

3=1 

and one may check that 

F(u) = Fi(u) + F 2 (oj). (6.16) 

The coefficients, which are needed in the calculation of 
F(ui), are list here 



hi 
hi 

fkl 

d n 

fl2 
fkl 



O21 - C12C22 sin uji 



1 — sin 2 u>i 

1 — C22 sin 2 uj\, 
02k 

1 — c 2 2 sin 2 uj\ 



-Ci 2 sinwi, 
-c n smui 2 , 



(6.17) 



\ 



1=1 

Qifc — /l*l/fcl 

i-ELiMhI 2 -I/hI 



with k=3, 4. 



C. an application of the POVMs 

In a recent work, Raynal et al. considered the ques- 
tion, which came from the implementation of the BB84 
by using the four quantum optical coherent states {| ± 
a >,| ± ia >} [35], of how to discriminate the following 



Pi = >< a| + I - a >< -a|), 

P2 = -{\ia >< ia\ + \ — ia >< —ia\), (6.18) 

and the authors expressed the optimal failure probability 
in terms of the mean photon number: 



F opt = e - M( 



COS^i| 



sin/z|), 



(6.19) 



with /i = |a| 2 . 

Here, we shall reconsider this problem with the 
POVMs in Eqs. (13-14). Writing the all the known pa- 
rameters in terms of a, we designed a program to get the 
optimal velue of F(u) in Eq.(6.15) by scanning in the 
parameters space < wi,o->2 _ ^ 27r. The final result of 
our calculation and the analytic solution in Eq.(6.19) are 
both presented in Fig. 9, while the optimal values of the 
u>\ and W2 are given in Fig. 10. Although in some regions 
of Fig. 9, small discrepancy still exists, our numerical cal- 
culations are consistent with the analytical solutions well 
in most parts of the parameter space, this fact shall give 
great supports to our present proposal. 



VII. DISCUSSION 

In the present paper, we always adopt a naive under- 
standing of the mixture: suppose a mixtures is denoted 
by p, for examples, p = Ylj=i >< ^ili m each run 

of the experiment, the input for our OPI is still a pure 
state belonging to the set {\^j >}. 

For the system of N=2 and N=3, we have given a se- 
ries of derivations to show why these cases can be solved 
by applying the argument of filtering, certainly, within 
the DTR. The case discussed in Sec. VI, is an important 
case of N=4. There are still others types of operations 
for N=4 system: (a) filtering | \Pi > from | \E f 2 > and 
1^3 > with > as the background; (b) filtering l^i > 
from |<I>2 > with > and |^4 > as the background; 
(c) discriminating three pure states with > as the 
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FIG. 11: The OPIs for the following cases: (a) filtering l^i > from > and |^3 > with |^4 > as the background; (b) 
filtering |$i > from |\]>2 > with |\]>3 > and |\1>4 > as the background; (c) discriminating three pure states with |\I>4 > as the 
background; (d) discriminating four pure states and (e) discriminating |^i >, |^2 > and {\^3 >, |^4 >}• 



background; (d) discriminating four pure states and (e) 
discriminating |*i >, |* 2 > and {|* 3 >, |\J> 4 >}. Their 
OPIs are given in Fig. 11. In principle, all these cases can 
be solved by reducing to filtering. For the cases where 
more than once filtering is needed, it's hard for us to find 
analytical solutions for the optical values. 



A important profit of our scheme should be mentioned 
here: the POVMs for each case shall be able to, although 
maybe not in a optimal way, complete the task of discrim- 
inating when the a priori probability of each state is not 
completely decided. 



In end of this paper, we would like to emphasize that 
: first, for a given case, if one could prepare the input in 
the one-photon state, then our OPI can be directly used 
for the optical experimental realization. Then, although 
the POVMs are from the one-photon picture, yet they 
are are general and state- type independent. Finally, a 
proposal, rather than a complete proof, has been given 
here in order to find a solution to the problem of the 
quantum state unambiguous discrimination. It's still a 
open question that: whether the task of quantum state 
unambiguous discriminating, either of pure states or of 
mixtures, can be solved by reducing it to the problem of 
quantum state filtering? 
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APPENDIX: SOME KNOWN MATRICES FOR 
THE LOW-DIMENSIONAL CASES. 

For N=2, 




with 

dct(0- 1 (2)) = 1 - |o 12 | 2 , h=t 2 = Vl-M 2 (1-3) 
For N=3, 

(1 O12 O13 \ 
021 1 23 , (1.4) 
031 32 1 J 

with A(3), the adjoint of 0(3), to be 

1 — |c>23 1 2 Ol3°32 — O12 O12O23 — Oi 3 \ 
O23O3I - 021 1 — 1 13 1 2 O21O13 - 23 , (1.5) 
O32O21 - O31 O31O12 - O32 1 — |oi2 | 2 / 
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The inverse 0_i(3) could be given by 0_i(3) = The overlap matrix , ± (3), is known with the form: 
i(3)/det(0) with 

det(0(3))= l-|o 12 | 2 -| 0l3 | 2 -| 023 | 2 

+012023031 + 021013032 (1.6) 

which has applications in 

/ dct(O) _ / dct(O) _ / det(O) 

h V r~h^F' <2 " V T~R^F '* 3 " V RTF' 

(1.7) 



013032 —012 



012023 — »13 



± (3) 



023031— 021 



V(l-|°23| 2 )(l-k 
032021— Q31 



V(l-|o 23 | 2 )(l-|oi3| 2 ) X /( 1 -|°23| 2 )(l-|oi2) 
^ 021013 — °23 



1 "12 —Q32 
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